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O . Abstract 

o 

Q>^ \ Influence of weak nonmagnetic impurities on the single-particle density of 

states p{uj) of two-dimensional electron systems with a conical spectrum is 
studied. We use a nonperturbative approach, based on replica trick with 
subsequent mapping of the effective action onto a one-dimensional model of 
interacting fermions, the latter being treated by Abelian and non-Abelian 



Q , bosonization methods. It is shown that, in a d-wave superconductor, the 

density of states, averaged over randomness, follows a nontrivial power-law 
behavior near the Fermi energy: p{ijo) ~ |w|". The exponent a > is calcu- 
^ ■ lated for several types of disorder. We demonstrate that the property /9(0) = 

is a direct consequence of a continuous symmetry of the effective fermionic 
model, whose breakdown is forbidden in two dimensions. As a counter exam- 
ple, we consider another model with a conical spectrum - a two-dimensional 
orbital antiferromagnet, where static disorder leads to a finite /9(0) due to 
breakdown of a discrete (particle-hole) symmetry. 
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I. INTRODUCTION 

The problem of weak static disorder in electronic systems with extended Fermi surface has 
been extensively studiedtl. It is well understood that randomness has a profound influence on 
the transport properties and weakly affects the density of states. The latter is no longer true, 
if the density of states (DOS) of a pure system vanishes linearly at zero energy: pioj) ~ Ic^l. 
In this case the standard procedure of averaging over impurities is complicated by the 
appearance of logarithmic singularities in the perturbative expansion of the single-electron 
self energynu. 

S(t^„) = iujn{-g\n{A/\ujn\) + •••)• 

Such a situation can be realized in two-dimensional systems with a Dirac-like, conical spec- 
trum E'^{p) = vIpI + VyPy, describing low-energy states near a degeneracy point (node), as 
well as in three-dimensional polar superconductors, where the order parameter has a nodal 
line. In both cases the pure systems exhibit a T^ low temperature specific heat. Since the 
exponents of thermodynamic quantities are used as an experimental criterium for selection 
of possible order parameters, it is important to study their stability with respect to disor- 
der. For 3D polar superconductors, the problem is simplified by the fact that, due to the 
finite size of the Fermi surface, the diagrams with crossing impurity lines are not divergent, 
and one can easily sum the remaining diagrammatic series. This was done by Gor'kov and 
KaluginH, who showed that weak nonmagnetic disorder gives rise to a special energy scale 
Qo ~ Aexp[— const /c], c being the impurity concentration, below which the DOS becomes 
finite p(0) ~ Qq. 

The renewed interest in quenched disorder in 2D systems^'^ has been stimulated by a 
number of experimental indications that the pairing state in the copper-oxide superconduc- 
tors may be of dj.2_y2 symmetry. In a recent publication, LeeQ applied the results of Ref. ^to 
analyze the role of nonmagnetic impurities in a 2D d-wave superconductor. He argued that, 
since a system with a nonzero DOS formally resembles a metal with a finite Fermi surface, 
the standard scaling argumentsll would indicate localization of all low-energy states. 

In this paper representing an extended version of our short reportQ, we show that, when 
averaging over impurities in 2D systems with a finite number of isolated nodes, the logarith- 
mic singularities appear in all self-energy diagrams, including the crossing ones. No specific 
subset of diagrams can be selected from the whole self-energy series expansion. What one 
finds here is a typical "logarithmic" situation when logarithmic singularities appear in all 
self-energy diagrams, including the crossing ones. Approximations based on partial sum- 
mation of the ("rainbow") diagramgHla do not apply. We overcome this difficulty by means 
of the field-theoretical approach based on replica trick. The static nature of the disorder 
allows to represent the effective, disorder-free action in terms of a one-dimensional model 
of interacting fermions. Using then bosonization technique, we demonstrate that, instead 
of creating a finite DOS at u; = 0, in the 2D d-wave superconductor the disorder changes 
the exponent of the DOS: p{ijj) ~ |ci;|", (0 < a < 1). The magnitude of a depends on the 
type of disorder. For a slowly varying (quasiclassical) random potential, when only forward 
scattering is present, 1 — a ~ c. The internode backscattering processes strongly reduce a 
at energies ~ Qo, making it concentration independent: a = 1/7. The picture indeed be- 
comes very close to the one with a finite DOS, but still not the same. In this situation, the 



renorinalization group for the conductivity can lead to an intermediate fixed point, whicli 
suggests a finite conductivity. 

Superconductors are not the only systems where a conical spectrum and, as a con- 
sequence, linear DOS can appear. Such a spectrum occurs in zero-gap (degenerate) 
semiconductorsB, and in hetero junctions where the contact is made between semiconduc- 
tors with inverted symmetry of bandgSo. It can be realized in 2D graphite sheetsc9, as well 
as for lattice electrons in a strong magnetic field, the most well known example being a 
tight-binding model of fermions with 1/2 of a magnetic flux quantum per plaquette (flux 
phase )lii. The conical spectrum is also a property of hypothetical orbital antiferromagnet 
and spin nematic statesllj. 

Despite the common feature of all these systems, - the existence of the degeneracy points 
(nodes) in their spectrum, the role of disorder can be quite different in each particular case, 
depending on the symmetry of the pure system. The result is determined by the structure 
of random terms which appear in the effective massless Dirac HamiltonianiS, describing low- 
energy states of the pure system. In the 2D d-wave superconductor, nonmagnetic impurities 
add to the Dirac model random gauge fieldd^, abelian or non-abelian, depending on the 
type of included scattering processes. The random gauge fields do not break particle-hole 
symmetry in the vicinity of the nodes and lead to a critical, power-law behavior of p{uj) at 
u; — >• 0. As we demonstrate below, this is a direct consequence of the fact that at u; — »• 
the equivalent ID fermionic model possesses continuous replica and chiral 7^ symmetries, 
remaining unbroken in two dimensions. 

On the other hand, in 2D systems with a particle-hole condensate, such as the flux phase 
or orbital antiferromagnet, the situation is different. The nonmagnetic disorder gives rise to 
a random chemical potential and random charge-density wave, the latter being equivalent 
to a random relativistic mass of the Dirac quasiparticles. We argue that the random Dirac 
mass alone would also lead to a critical behavior with p(0) = 0. However, the random 
chemical potential, which is always present even for a weak quasiclassical disorder, breaks 
the particle-hole symmetry and drives the system away from criticality. The special role 
of the random chemical potential is that it lowers the continuous chiral symmetry of the 
corresponding ID Fermi model down to a discrete one. The spontaneous breakdown of the 
latter is signalled by a finite p(0), as previously pointed out by Fisher and Fradkinllll. It 
is worth noticing that, for a d-wave superconductor, the appearance of a random chemical 
potential in the effective Dirac model could be only caused by random magnetic fields. 

The paper is organized as follows. In Section II we formulate a model of two-dimensional 
d-wave superconductor with non-magnetic impurities. In Section III we give a perturbative 
analysis of the disorder in terms of the diagram expansion for the single electron Green's 
function. The expansion demonstrates that the number of singular diagrams in the given 
n— th order grows with n, and the problem is non-perturbative. In Section IV we average 
over the randomness using the replica trick and reformulate the replicated model as a two- 
dimensional Euclidean relativistic field theory. In Section V we consider a simplified version 
of this theory, namely, we consider only scattering within each node. This simplified theory 
is solved by the ordinary bosonization. In Section VI we consider the scattering between 
nodes; the corresponding model is solved by the non-Abelian bosonization. In Section VII 
we describe another model with a conical spectrum - a 2D orbital antiferromagnet. We 
claim that this model has a finite p(0). The paper has an Appendix where we discuss some 



properties of the Bethe-ansatz solution of the rephcated model of impure d-wave supercon- 
ductors. 



II. A MODEL OF IMPURE D-WAVE SUPERCONDUCTOR 

We start from the following model describing low energy electronic degrees of freedom 
in a 2D d-wave superconductor: 

H = E<k)c+,Ck. - E[^(k)4T4i + H.c.]. (1) 

k,CT k 

Here e(k) is a tight-binding spectrum in the normal state which respects full point symmetry 
of the underlying (square for the copper oxides) lattice, and A(k) is a pairing amplitude, 
odd under 7r/2-rotations in the k-space (the d^2_y2 symmetry). Without loss of generality 
we can choose e(k) = —2t{coskx + cos ky) — /i and A(k) = Ao(cosfcj; — cos ky). It will be 



assumed that Aq ^ t. The quasiparticle spectrum E(k) = ±-y/e^(k) + A2(k) has four nodes 
at the Fermi level which we denote 1,1,2,2 (see Fig. [l|). Their positions are ki = — kj = 
(fco,/co); k2 = —kg = (— fcoj^o); where k^ = arccos(/i/4t). Associated with these nodes are 
gapless excitations with formally relativistic (conical) spectrum at |a;| -C Aq. The single- 
electron density of states at these energies goes linearly with u: p{uj) ~ |ct;|, giving rise to 
extra T-factors in thermodynamic quantities of pure d-wave superconductors. 
Let us introduce the Nambu spinor 

The Hamiltonian (1) in this notations is given by 

H = Y. "^k [(^k - ^An + A(k)r2]$k, (2) 

k 

where Tj are the Pauli matrices. We linearize the spectrum close to the nodes and pass to 
a continuum description in terms of four Fermi fields ipj{x), j = 1, 1, 2, 2. We choose a new 
coordinate system rotated by angle n/A with respect to the original one. It is also convenient 
to make a tt/2 rotation in the Nambu space about the r2— axis. The resulting continuum 
model is then given by the sum of four Dirac-like Hamiltonians: 

H= E fd^xtlj+{x)Hj{x)^/jj{x), (3) 

j=l,l,2,2 

Hi{x) = -Hi{x) = -ividiTi - iv2d2T2, (4) 

H2{x) = -H2{x) = -iVid2Ti - iV2diT2, (5) 

with ^1 = 2\/2ta sin kQ,V2 = "\/2Ao sin ko- 

Processes of scattering on nonmagnetic impurities are described by 

Himp = J2 d'^xVjj>{x)iljj^{x)Tiiljf{x), (6) 

jj' 



where Vjj'{x) = V*,j{x) are random fields with Gaussian distributions. Since Vjj'{x) are 
wave packets with average wave vectors kj — k^v, not all of them are independent. For the 
square geometry presented on Fig. [1], we have the following independent components: 

(i) the real field Vjj{x) = Vo{x) representing a slowly varying (quasiclassical) component 
of the impurity potential; 

(ii) the complex fields Kii(x) = Vi(x) and 1^22(2;) = V2{x) representing a back scattering 
with the momentum transfer 2ko', 

(iii) the complex fields Vuix) = V^i = W\{x) and Vyii^x) = V2i(x) = W2{x) corresponding 
to transitions between nearest nodes. 

In what follows we assume that all these fields are distributed according to the Gauss' 
law: 

<Voix)Voix')> =\o5^^\x-x'), 
<V,ix)V;ix')> =X,5^^\x-x'), (j = l,2), 
< Wj{x)W*{x') > = X26'^^\x - x'), (j = 1,2). (7) 



III. SELF ENERGY DIAGRAMS 

For a given realization of static random potential, the one-particle Green's function 
tensor written in the mixed x — Un representation {un = (2n + l)7r//5 being the Matsubara 
frequency) satisfies the following integral equation 

Gjj'{x, x'; Un) = Sjf G^j{x - x; Un) 

+ X! / '^^^" ^^(^ ~ ^"' ^n)Vjjii{x") fi Gj»j'{x", x'; Un), (8) 

where 

G'^jix - x'; ujn) = [luJn - Hj{x)]-^6{x - x') (9) 

is the Green's function for the pure system. 

For weak electron-impurity scattering, the standard procedure (see, e.g. 0) consists in 
developing a perturbation theory expansion in Eq.(|]), with subsequent averaging of products 
of the Gaussian random fields Vjj'{x), appearing in each term of this series. The averaged 
Green's function 



< Gjji{x, x'\ ijJn) >= 6jj>Gj{x — x'; Un), (10) 

written in momentum representation, satisfies the Dyson equation 

Gj^{k,LJn) = iuJn - Hj{k) - f:j{k,Un). (11) 

Here 

^i(k) = -Hi{k) = kmn + k^v^h, (12) 



^2(k) = -i^2(k) = k2Vin + hV2f2. (13) 

Sj(k, C(j„) is the self-energy operator incorporating all effects of the electron-impurity scat- 
tering. In what follows, we shall analyze lowest-order diagrams for Si. In calculations, a 
cutoff prescription, defined by conditions |ifj(k)| < Aq, will be used. 

In the lowest (Born) approximation Si is given by the first diagram shown in Fig. ^j. It 
is momentum independent and contains a logarithmic singularity, as previously pointed out 
in a number of papersoQ: 

+A2(G^(k',^„)+G'^(k',,;„))]fi 

d^p 1 



-iuJn {93 + 91 + '2 92) / 



p|<Ao(27r)2p2+^2 

. 93 + 91 + 292 Ao , . 

-i^n 7, In - — -. (14) 

ZTT \U!„\ 



Here 



Ai,2 Ao , . 

91,2 = , ^3 = (15) 

fit;2 t;it;2 

are small dimensionless coupling constants, proportional to the impurity concentration c. 

In fact, logarithmic singularities appear in all orders of perturbation theory. Consider 
second-order diagrams shown in Fig. ^j. The "rainbow" diagram is obtained from the Born 
one by first-order renormalization of the internal electron line. The node indices ji and j2 
take arbitrary values from the set (1, 1, 2, 2). The result is easily found to be 



V(2a)/ N • (93 + 91+^92) r , 2 -^0 , ^n 
S^ >{uJn) = -lUJn 777-T^ In -. r + 0(1 

(27r)2 \ujn\ 



(16) 



The diagram with crossing impurity lines, shown in Fig. ^, corresponds to a vertex 
renormalization in the Born self-energy diagram. The momentum conservation restricts 
the allowed distribution of the node indices {ji,J2,J3) by the following nine possibilities: 
(111), (111), (111), (221), (122), (221), (122), (2l2), (2l2). Using power counting and symme- 
try arguments, one can show that, in each diagram, the maximum power of log |a;„| appears 
with the prefactor ~ iuJn- This allows to simplify calculations by setting the external mo- 
mentum k = 0, in which case 



S(2^)(k = 0,^„)|(,,,,,3) =/ 



£k (fq 



(27r)2 (27r)2 
xfi ^"^(k,^^) fi G';^(q,^„) fi G;3(q-k,^„) fi. (17) 

Consider, for example, the impurity scattering in the vicinity of a single node, corresponding 
to the case (111). Here we are able to rescale all the momenta to get a circular cutoff in 
each integral: 
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^'?U-n) = -^-n9l J J--J 



(Pk f (Pq 



(|fc|<Ao) (27r)2 J(\q\<Ao) {27t] 



2 



k.(q-k) 

Integrating over the angle between k and q yields — 27r0^^ ijp — q^), where 0a;(x) is a smeared 
step function, with the width of the order |u;| at s = 0. The remaining integration over the 
region |k| > |q| gives a log \uJn\ contribution: 

(26) . gl f^o qdq f^o k dk 

VDl'^nj-^^n ^277)2 io ivl + q^J, ul + k^ 



lUJr, 



9i ,2 A 



2(27r)2 \uj, 



In' Pt- (19) 



Quite similarly 



fi'35'l , 2 Ac 



^ iii)(^") + ^ iii)(^") - -'^^n jTT-T^ In 1 r- (20) 



In cases when intermediate states belong to neighboring nodes, as it takes place for di- 
agrams (221), (122), (221), (122), (212), (212) the situation is somewhat more complicated 
due to the difference between the two velocities, vi ^ V2. The two ellipses, representing 
the surfaces of constant energy for states 1 (1) and 2 (2), are strongly elongated in per- 
pendicular directions and, therefore, cannot be simultaneously transformed to circles by a 
global rescaling of the momenta. The available phase space is then reduced by the amount 
determined by the anisotropy parameter 

7 = t^sM ~ Ao/t < 1. (21) 

As a result, an additional energy scale 7A0 enters the problem, and the corresponding 
diagrams become dependent on the ratio |ci;„|/7Ao. For example, 

(26) / X , v(26) r, \ - o.-, , „2 f d'^k f (Pq 



(212) v"^"/ + ^("212) v"^"/ ~ 2zCo'„ (^2 



(2i2)V-n; -r -(212) V-n; " -^-n iJ2 7„„^^„, /o^^2 /„„,,_^„ ,„„,^^„, (0^\2 



(|A:|<Ao) (27r)2 7{|gi|<7Ao,k2|<Ao) (27r) 
(a;2+F)(^2+^-2^2+^2^2^[^2 + (q_k)2]- I ) 



Integrating first over k, we get 



^(26) ( . , ^(26) ( . . '^iujn 2 H^^ Z"^" > 1^ [Ao/maa;(|^„| , |g|)] 

4ii)(-n) + S^(-n)--4^^^2X d^^i dl- ,2 +^-2^2 ^^2^2 • (23) 

At energies 7A0 < \ujn\ <^ Aq, the integral in |^ is mostly contributed by regions gi -C 
7A0, 7A0 <^ g2 ^ Ao, and is nearly independent of Un'- 



.(26) ^, , ^ , ^(26) I,, ■, ^ ^^ 

{2^y 



Egy^„) + Sgi)(^n) ^ -7^ 7A0 [1 + 0(7 In 7)] sign(.;„). (24) 
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The log^-term only appears at lower energies, |ti;„| ^ 7A0, in which case integration over 
|g| < 7A0 yields 



.(26) / N , ^(2b) / N . 5'2 1 2 7^0 



(^„) + ^fU^n) ^ -luJn 77?^ In^ f^. (25) 



For the remaining four diagrams ((221), (122), (221), (122)) the effect of the anisotropy 7 is 
even stronger. It can be shown that at |a;„| -C 7A0 these diagrams contain only first power 
of the logarithm: S ~ zco'„(735f2 ln(7Ao/|ct;„|). 

The above analysis shows that, in 2D electron systems with a conical spectrum, impurity 
scattering drastically differs from that in a normal metallic state with a finite Fermi surface. 
The structure of the perturbation theory expansion for the self energy reveals a "logarithmic 
situation" resembling mass renormalization in models of interacting fermions in one space 
dimension. (This analogy, being realized in earlier paperalllla, will become more transparent 
in subsequent sections, where the effective field theory in (1 + 1) dimensions is discussed). 
Notice that, after analytic continuation, iujn — > ciJ + i5, the leading logarithmic singularities 
occur in the real part of S(a;). This means that, in systems with a conical spectrum, 
disorder mostly affects renormalization of the single-particle spectrum, as opposed to the 
usual picture in systems with a finite DOS, where life-time effects dominate. 

Actually, the logarithmic singularities appear in all orders of perturbation theory, includ- 
ing crossing diagrams. In two- and higher-dimensional metallic systems with a finite Fermi 
surface, such diagrams are known to be relatively smallli^l. In the case of weak disorder, 
kpl ^ 1, Z being the mean- free path, strong restrictions over the relative angles between the 
momenta of intermediate one-electron states effectively reduce the available phase space, 
making diagrams with crossing impurity lines proportional to powers of the factor (l/kpl)- 
However, in ID systems, where the Fermi surface is represented by two points, ifej?, no such 
reduction is possible, and all diagrams are equally important, as is known from the theory 
of one-dimensional localizationoEj. In this respect, the electron- impurity scattering in 2D 
systems with conical spectra resembles the ID case: the crossing diagrams also contain log- 
arithmic singularities and therefore should be treated on the same footing as the rainbow 
ones. For a finite number of degeneracy points (nodes), no selection of diagrams is possible. 
An approximation, in which one confines consideration to the class of rainbow diagrams 
and then solves self-consistently the resulting Dyson equation for S(u;)eI, can be justified 
only in the limit of (infinitely) large number of nodesB, or when one extends consideration 
to a higher dimensional case, e.g. 3D polar superconductors with a line of nodescl. This 
approximation leads to a finite DOS at a; = 0. 

In the next sections, all singular self-energy diagrams will be taken into account, using 
replica trick. Treating the effective, disorder-free field-theoretical model by abelian and non- 
abelian bosonization, we shall show that, for a weakly disordered 2D d-wave superconductor, 
the low-energy DOS follows a power-law behavior. The corresponding critical exponent 
depends on the type of electron-impurity processes included into consideration. 



IV. REPLICA TRICK AND EFFECTIVE TWO-DIMENSIONAL FIELD THEORY 

In what follows, we shall be mostly concerned with the forward and backward scattering 
processes, putting A2 = 0. In this case the nodes (1, 1) and (2, 2) can be considered indepen- 
dently. Without loss of generality we can choose the (1, 1) pair. Rescaling the coordinates 
and the fields 

x^fixi, y ^ V2X2, 

iji{x) -^ — v^i(x), v^i(x) -^ — iJi{x) 

^VxV2 ^VxV2 

we arrive at following generating functional for the (unaveraged) Green's function with a 
fixed Matsubara frequency ujn- 



Z[A] = / Df]D'rjexp{-S[A]), (26) 

S = (fxf]{x)[{-idi + A3(x)(T3)ri 

+ {-id2 + Ai{x)ai + A2{x)a2)T2 - icu„]r/(x). (27) 

Action ( p7D describes two-dimensional massless fermions interacting with a random non- 
Abelian gauge potential. The Pauli matrices ctj act on the isotopical indices (1, 1), and the 
probability distribution of the gauge field is 

P[A] = I DA exp[- fd^x^^]. (28) 

The matrices T2 = 70 and ti = 71 form a representation of the two dimensional Clifford 
algebra, i.e they are the Dirac matrices of our problem. The new fermionic fields are related 
to the initial ones as follows: 

f] = —iip^i 7] = ip. (29) 

In terms of the new fermionic variables, the DOS is given by 

p{uj) = Re[Tr < r/(x)r/(x) >]\i^„^^+i5. (30) 

TTViV2 

We average over the disorder using the standard replica trick. For this purpose we 
replicate the action ( [27| ) r-times (r is the number of replicas) and integrate the fermionic 
partition function with the replicated action over the gauge fields. The result is 

Zr = Df]Driexp{—Sr), (31) 

+ l^[{f]a,all(rlpVa,pf + (^a,a7l^a/3„,^)2}, (32) 



where the Greek indices are isotopic and the Latin ones are reserved for rephcas (a = 
1, ...r). It will be convenient to rewrite the problem in the Hamiltonian formalism treating 
Xi asimaginary time and X2 as space coordinate. The quantization rules are f/70 = t]^, where 
the spinors {rjR, tjl) and their Hermitian conjugate satisfy the standard anticommutation 
relations. The Hamiltonian is given by 



H = dx{r]^^^{-id:^Ts + u;„r2)?7a,, 



-y(: Vta^3(^lpVa,f3 -Y + -^ [{■■ Vta(^il3Va,/3 -f + -j{- Va,a(^lpVa,f3 -f] }, (33) 

where the Matsubara frequency u;„ plays the role of relativistic mass. 

Notice that the Hamiltonian (|33| ) possesses continuous SU{r) replica symmetry. At a;„ = 
it is also invariant under chiral rotations of Fermi fields: (the continuous 7^— symmetry) 



T] —^ exp{iv^(p) rj. As follows from (|30D, a finite p(0) would mean the existence of a nonzero 
order parameter < fjri >, indicating a spontaneous breakdown of these symmetries. However, 
continuous symmetries cannot be broken in two dimensions; therefore p{uj) should vanish at 
cu — s> 0. This is confirmed by direct calculations presented below. 



V. MODEL WITHOUT BACKSCATTERING. BOSONIZATION 

In this Section we shall focus on the slowly varying (quasiclassical) part of the random 
potential, i.e. we set gi = §2 = 0. In this case the problem is reduced to the single node 
one. The isotopic indices can be suppressed, and the resulting model can be treated by the 
standard abelian bosonization methods . 

The bosonization rules are well knownO: 

'jR,a + <^L,a = —J= dx4>a , JR,a — JL,a = 1= Ha; 

'IT \/TT 



■ VR,aVL,a + vtaVR,a ■ = "A COs(V47r0„). (34) 

Here 

jRiL),a =■■ VR(L),aVRiL),a ■ (35) 

are the current operators for right (left )-moving particles with the replica index a; (/)a(a;) and 
H(i(x) are scalar fields and their conjugate momenta, respectively, satisfying the canonical 
commutation relations: [<^a(a;), nf,(y)] = i6abS{x — y); A ~ Aq is the ultraviolet cut-off. The 
Bose model equivalent to (BBf) reads: 



H= fdx{W[^l + {d,^a)']-^J2^anb-UJnAj2cOsiV^<Pa)}- (36) 

In order to proceed further we need to diagonalize the quadratic form J2ab ^a^b = 
H^MH, where Mab = 1 for all a and b. We do it with an orthogonal transformation U: 
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u = up, u^u = r, 



and obtain 



Uar = ^ (a = l,...r). (37) 



H = Jdx{\ E [P'a + (9.0.)^] + ^ [(1 - "-z^)P^ + {dM'] 

a=l 
-LUnAj2'^Os{V^Uab(l)b)}- (38) 



The last step consists in rescaling 0^, and Pr-: 



(f)r — > 70r , Pr -^ -Pr, 

7 
^2^(l_^)i/2. (39) 

TV 

The resulting theory is described by the Hamiltonian 

/ '' y 

d^iH -flPa + {d.<Pa?] - UJnA E COs(v^K60fe)}, (40) 

a=l ^ a 

where 

Ua = l ia< r), Ur =7, 

1 

Vab = Uab (a, & < r), Kr = -• 

r 



The critical dimension of the cosine term in Eq. (^Oj) is given by 

Aa = EKl = l + ^^^- (41) 

We see that it does not depend on a and remains well defined in the replica limit: 

A ^ A,(r ^ 0) = 1 - 1^. (42) 

ZTT 

Simple scaling arguments allow us to estimate the DOS at finite u. First we note that 



in Eq. f^OD the cosine term is a relevant perturbation generating a new energy scale a; (or, 



equivalently, correlation length ^c ), given by 

l-A 

Cj r^C^ ^uj\--\ . (43) 

\M) 

This formula correctly reproduces perturbation theory expansion for the self-energy operator 
in the single-node approximation [gi = (?2 = 0): 
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^jtUJn) I _ 1 _ ^ _ 1 (A\ "^ 

I ja;,i— »a;+j(5 ^ I i i I 

lUJn UJ \\^\J 

27r |cj| 2 \2txJ \oj\ ^ ' 

in agreement with the estimation of first- and second-order diagrams (see Eqs.(P^,(P^ 
and (|19D). Applying then bosonization rules (Q) to Eq.(|30D we establish that piu) ~ 
{cos{'\f^Vab4>b)) ■ This average is proportial to i~^. Using then (^3[), we find a power-law 
behavior of the low-energy DOS 

p{uj) ~ |tu|" (45) 

with a nonuniversal exponent 

A 1 - ^3/27r 

a = r = 1 — < 1, (46) 

2-A l+^3/2vr ^ ^ 

depending on the impurity concentration. 

VI. IMPURITY SCATTERING IN MODEL WITH SEVERAL FERMI POINTS 

In this Section we consider the model where scattering processes transfer electrons be- 
tween the opposite Fermi points. The corresponding replicated model was described in 
Section IV, Eq. (|33|); here we rewrite it explicitly in terms of chiral fermions: 

H = Ho + Hi + H2, (47) 

Ho = J dxT]^^^{-id^T3 + UJnT2)r]a,a, (48) 

H^ = fdx{-2gsi: Jpl : + : Jpl :) + 2g, ^ (: X^X^ : + : jyi :)}, (49) 

•' i=l,2 

H2 = AJdx[g, : Jpl : +g, Y. ■ Vl :], (50) 



4 = 1,2 



where 

ji + ^aP 

-Jr - VR,aa—VR,al3, 

Jl = Vtaa^VL,aP (51) 

are chiral components of the SU{2) currents. 

We shall apply to the model (J^7\ ) the procedure of non-Abelian bosonization developed 
by WittenEj. The approach is based on the following key moments. The first one is that 
the Hamiltonian of free fermions ( ^HD can be expressed in terms of current operators (the 
Sugawara constructionllZHl3) . The Sugawara form of the Hamiltonian of free massless fermions 
with a general U{1) x SU{N) x S'[/(r)-symmetry (in Eq.(|D N = 2) is 

12 



= / dx[—{: Jr{x)Jr{x) : + : Jl{x)Jl{x) :) 

i\ + r ^^^ 

E(: Jr{^)J%{x) : + : J£(x)JE(x) :)], (52) 



N + r 



a=l 



where J, J* and J*^ are the f/(l), the SU{N) and the SU{r) currents defined as in Eg. (pl|) , 
but with generators of the corresponding algebras instead of the Pauh matrices. Gn and 
Gr are the total number of generators of the su{N) and the su{r) Lie algebras. The second 
key moment is that the currents from different algebras commute. Therefore Hq is a sum of 
three mutually commuting operators. 

Now notice that the interaction terms (^) , (pOD contain only spin currents and therefore 
do not affect the spectra of 5'[/(2)-singlets. Since this interaction is attractive, the spectrum 
in the SU{2) channel becomes gapful and, as we show later, this gap persists in the replica 
limit. Therefore, if ujn = 0, the spectrum well below the gap Qq is described by the rest of 
the Hamiltonian {^7\), in other words by -f^e// = H^,^^ + H^^,y The Hamiltonian H^^ 



IS 



the Hamiltonian of the Wess-Zumino-Witten model. Its spectrum is the subsector of the free 
fermionic spectrum generated by the SU{r) current operators. The model is conformally 
invariant and exactly solvablecj. Its primary fields transform according to representations 
of the SU{r) group. In particular, the field from the fundamental representation, which we 
denote gab, has the following scaling dimensioned 

A. = ^- (53) 

We see that A^ is singular at r — > 0, which indicates that the (/-field cannot appear alone 
in the physical sector. In fact, as we shall see, a local bilinear form of two Fermi operators, 
f]ri^ which appears in the "mass" (i.e. cUn-proportional) term of the Hamiltonian, as well as 
in the definition (30) of DOS, is represented by a combination of g and a phase exponential 
of the scalar (t/(l)-symmetric) field 



/47r 
riT] = Q r^ Qo Re [^exp(iW— 0)]. (54) 

The scaling dimension of the exponential is equal to 1/Nr. So, the two singularities cancel, 
and the resulting scaling dimension of the operator Q is finite in the replica limit: 

1 r — 1 /r 1 

An = lim [— + ^1 = ^T. (55) 

It is instructive to repeat the above derivation in a more traditional form introducing the 
replica tensor Qab- For the sake of simplicity we assume gi = gz- Let us rewrite the action 
(B3) in the form suitable for the Hubbard- Stratonovich transformation: 
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~g[{Va,aVb,a){Vb,l3Va,p) - {Va,al5Vb,a){f]b,f3j5Va,l3)] + Sirr, (56) 

where 75 = iri and Sirr contains terms with interactions near the same Fermi point. Such 
terms renormahze velocities and are not important for the current discussion. Now we intro- 
duce the auxilary field Qap and decouple the interaction term by the Hubbard- Stratonovich 
transformation. Integrating formally over the fermions, we obtain the partition function for 
the replica tensor field: 

Z= f DQ+DQexp{- j <fxL), 

L = ^Tr{Q+ - ujI){Q - ul) 

-N Tr\n[i^^d^ + (1 + 275)^/2 + (1 - ^75)Q+/2]. (57) 

Following the conventional wisdom, we shall look for a saddle point of the exponent and 
then derive an effective action for fiuctuations around this saddle point. We assume that 
the saddle point configuration of Q is coordinate independent and as such can be chosen 
as a diagonal real matrix: Q{x) = diag{\i, ...\r). The density of effective action on this 
configuration is equal to 

Seff = V[i^^^i^ + ^ In ^]. (58) 

The saddle point value of Q, being a point of minimum of this function, satisfies the equation: 

Xn — UJ NXa ^ \Xa\ „ ,, ^, 

— + ^ln^^ = 0. 59 

g TT A 

At cu = we have the following solution: 

Xa = Aexp[-n/Ng] = Qo. (60) 

The existence of the saddle point does not mean, however, that the system has a non-zero 
order parameter (Q). Such order parameter in the present case is the density of states at 
the Fermi energy: 

dZ 1 

p(0) = limr-i — U=o = lim — Tr{Q + Q+). (61) 

1*^0 OUJ r~*Q gr 

The appearence of the average (Q) breaks the continuous symmetry U{r), and it is well 
known that such symmetry breaking cannot occur in two dimensional systems due to strong 
transverse fiuctuations of Q. A finite Re (Q) arises only at finite u when the transverse 
fiuctuations acquire a finite correlation length ^^ ~ uj^^^^'^^^\ A is the scaling dimension of 
the operator ReQ. In this case we have: 

Re {Q) r^ Qo uj^ . (62) 
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Now we are going to demonstrate that the described quahtative picture survives the 
rephca hniit and maintains a contact with the previous treatment. In order to show that 
the radial fluctuations of Q remain gapful even at r = 0, we use the exact solution of the 
model (^). This exact solution was found by one of the authorst^. It was shown that the 
model ( p6[ ) at cu = is a relativistic limit of the model directly solvable by the Bethe ansatz: 

H = j dx[i,X^{-\dl - ^4)^a,« - 9 (^:„^^«)(^.>a,/3)]. (63) 

The basic features of the exact solution coincide with those derived from the non-Abelian 
bosonization. In particular, it was established in Ref. |l^ that, in the relativistic limit, the 
spectrum of model ( ]63| ) consists of three sectors. Excitations of one sector are f/(r)-singlets 
and have a spectral gap given by Eg . (|60D . The low energy sector contains gapless excitations 
which are S'[/(A^)-singlets. One branch of gapless excitations consists of the U(l) scalar field 
described by the effective action 

S^ = \jd?x{d,ct>f. (64) 

while the other is described by the Wess-Zumino-Witten model at the critical point, i.e. by 
^su(r) ^^ Eg. (p2D . The latter can be also written in the Lagrangian form with the following 
action: 

f N N r'^ 

S2 = J d\{—Tr{d,g+d^g) + ^^ J^ d^e,,, Tr{g+dagg^d,gg+d,g)], (65) 

where g is & matrix from the SU{r) group {g^g = I,detg = 1). The second term in the 
r.h.s. of Eg.(^) is the so-called Wess-Zumino term. It has a topological origin. Despite of 
the fact that it is written as an integral including the additional dimension, its actual value 
(modulus 2mN) depends on the boundary values of g{x,^ = 0) = g{x) {g{x,^ = 1) = 0). 
This property follows from the fact that the Wess-Zumino term is proportional to the integral 
of the Jacobian of transformation from the three dimensional euclidean coordinates to the 
group coordinates, thus being a total derivitive in disguise. Excitations of the Wess-Zumino- 
Witten model are gapless and the correlation functions are well studied (see Ref. ^0]). The 
presence of the Wess-Zumino term is absolutely crucial for the criticality - the model without 
it undergoes a strong renormalization. 

Let us establish a connection between the exact solution and the semiclassical approach. 
The massive sector obviously corresponds to radial fluctuations of the order parameter field 
Q. The rest of the action (except of the Wess-Zumino term) represents the first term in the 
gradient expansion of the Tr In-term in Eg . (|57D . Indeed, at small momenta \p\ <^ Qq this 
term is equal to 



N 



j (fxTr{d^Q+d^Q). (66) 



%-kQ 

Now let us write down the Q- field in the form (Q). Substituting this expression into Eg.(|66D 
and taking into account that Trg^d^g = as a trace of an element of the algebra, we 
reproduce Eg.(p3) and the first term of Eg.(K3|). It comes not entirely unexpected that the 
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naive gradient expansion misses the important Wess-Zumino term: the latter is a Berry 
phase and requires special care. In order to avoid possible calculational difficulties, we have 
resorted to the exact solution. In this way we obtain the result having an independent value 



(see also Ref. 21 



-Tr Hi^^d^ + (1 + 275)^^0/2 + (1 - «75)^+Qo/2] = 

j d\{^Tr{d,g^d,g) + ^ f^ d^ e.,e Tr{g+dagg^d,gg+d,g)} + 0{Qf), (67) 

where g belongs to f/(r) group (r is now an arbitrary number not related to the number of 
replicas) . 

The fact that radial and transverse fluctuations are decoupled at the level of the Bethe 
ansatz equations allows one to consider the corresponding replica limits separately. The 
proof of the fact that the replica limit is well defined for the gapful sector is given in the 
Appendix. Here we consider only the replica limit for the correlation functions of Q. For 
the pair correlation function we have 



< Q{x)Q^{y) > =Ql< exp[i^47r/A^r0(x)]exp[-i^47r/Arr0(i/)] >< g{x)g+{y) > 

^ la; _ ^|-2(l/iV.+A,) ^ 1^ _ y\-(^/^'\ (68) 

in agreement with the above estimation of the scaling dimension of the operator Q (see 
Eg . (pSD ) . According to Eq. (|6^) the dimension 1/N'^ gives the following estimate for the 
density of states at |cc;| <C Qo^ 

p(cu)~cui/(2^'-i). (69) 

In fact, even in the case N = 2 the corresponding exponent is extremely small: 1/7 and the 
specific heat at low temperatures goes almost linearly: c^ ~ T^/^. 

So far the impurity scattering between the nearest-neighbor nodes has been neglected. 
When these processes are taken into account, the effective replicated action can be easily 
shown to have the form: 

Seff = Si + S2 + Si2, (70) 



Si= I d^x{r]iaa{VidiTi + V2d2T2 + UJn)r]laa 

1 1 

+ 253(^laanO-^/3^1a/3)^ + ^9l[{Vaia'r2(ripVlaf3f + (^aia^2C^a/3^1a/3)^]}, (71) 



5*2 = / d'^x{riiaa{Vld2Ti + V2diT2 + UJn)r]lo 



+ ^93{V2aarialf^r]2apy + -gi[{ri2aar2(ylpri2apf + (^2aar2a^/3r/2a/3)^]}, (72) 

5'l2 = 1^92'^ j d'^x{{f]iaaTl(ylprijalif - {Viaanrjjaaf 

+ {ViaaT2(^lf3Vjapf + {ViaaT2(^ipVjal3f}, (73) 
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where fermionic fields r]i and ri2 describe two pairs of opposite nodes, (1,1) and (2,2), 
respectively. The action (|70D is still symmetric under continuous 7^— transformations 

Vja ^ e^"^^ r]ja, ma ^ ^.a e*""^ (74) 

Therefore the DOS vanishes as u; — > 0, even when the nearest-neighbor internode scattering 
is included into consideration. This conclusion does not depend on the ratio Vi/v2. However, 
the crossover from the two-node critical regime to the four-node one will be shifted towards 
lower energies on increasing fi/f2. 



VII. NONMAGNETIC IMPURITIES IN 2D ORBITAL ANTIFERROMAGNET 

In this section, we consider effects of nonmagnetic disorder in another 2D electron system 
with conical energy spectrum - the orbital antiferromagnet state (OAF). This state was 
earlier suggested as one of possible phases of weakly interacting fermions on a square lattice 
with a half-filled energy bancO. The OAF phase can be viewed as resulting from anomalous 
particle- hole pairing < cj^^Ck+Q,a > with momentum transfer equal to the nesting vector 
Q = (7r,7r). Such a pairing breaks translational, point and time reversal symmetries, but 
preserves spin-rotational invariance. At the microscopical level, the OAF is characterized by 
nonzero local charge currents circulating around plaquettes in an alternating way. It actually 
represents a weak-coupling analog of the flux phase studied by Affleck and Marstoncj in the 
large-?/ Hubbard model. 

Since the OAF is a state with a spin-singlet particle-hole pairing, and since we wish to 
discuss effects of nonmagnetic impurities, spin degrees of freedom are irrelevant and will be 
suppressed in what follows. Then the mean-field Hamiltonian describing a 2D OAF state 
can be written as 

H, = E[e(k)c+Ck - zA(k)c+Ck+Q], (75) 

k 

where the tight-binding spectrum e(k) and the gap function A(k) have the same two-cosine 
structures as those for the d-wave superconductor (see (P) with /x = 0. The quasiparticle 



spectrum E{k.) = ±Je'^(k) + A2(k) has zeros at points ki = — kj = (vr/2, vr/2), k2 = — kj = 
(— 7r/2, 7r/2). However, due to unit cell doubling taking place in the OAF state, the reduced 
Brillouin zone corresponding to the new translational symmetry, coincides with the area 
enclosed by the original square Fermi surface ABCD, shown in Fig. ^ Therefore, there are 
actually two unequivalent nodes 1 = 1 and 2 = 2 in the OAF phase, with two associated 
branches of low-energy excitations with conical spectra. 

Assuming that q lies inside ABCD, we introduce a 2-spinor 

and rewrite Hamiltonian (75) in a matrix form 

Ho = T.K i<q)^3 + A(q)f2] %. (76) 
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Scattering on nonmagnetic impurities is decribed by 
Himp = ^V{k- k')CkCk' 

k,k' 

= Y. [K.(q - q'X*q' + ^i(q - q')*q^i*q'] (77) 



q,q 



where ^(k — k') is the Fourier transform of a Gaussian random potential defined on the 
original lattice. In the OAF state, this potential gives rise to two random fields, Vo(q) = ^(q) 
and ^i(q) = ^(q + Q), whose small- momentum parts describe the chemical potential and 
charge- density- wave (CDW) fluctuations, respectively. 

The physical difference between Vq and Vi allows to expect qualitatively different effects 
of these two types of disorder on the low-energy DOS in the OAF phase. Consider a local, 
slowly varying fluctuation of the chemical potential, such that Vo(x) is nearly constant in a 
macroscopically large region. The finite value of Vq will then (locally) break the particle-hole 
symmetry of the OAF state and open a finite Fermi surface, thus resulting in a finite DOS in 
the given region: Nioc ~ |Vo|. The latter conclusion does not depend on the magnitude and 
sign of Vo(x). Therefore, the DOS averaged over realizations of the field Vo(x) is expected to 
be finite at zero energy. On the other hand, local CDW fluctuations tend to open a gap in 
the excitation spectrum, so that, for V2(x) ~ const, the local DOS will vanish in the whole 
energy range \E\ < IV2I. Therefore, we expect in this case the averaged DOS to vanish at 
E -^ even faster than in the pure system. These expectations will be confirmed by the 
analysis presented below. 

Following the same strategy as in Sec. 4, we pass to the continuum limit based on sep- 
arating states near the nodes. For simplicity, we shall ignore scattering between nodes 1 
and 2. We shall argue later that the internode scattering does not alter main conclusions 
obtained in this section. Specializing to the node 1, we arrive at the Euclidean action 



S[Vi\ = d xri{x)[-idiTi - id2T2 + /i(x) - m(x)r3 - iujn]r]{x) (78) 

which describes two-dimensional fermions with a random mass m{x) = Vi{x) and chemical 
potential /i(x) = Vo(x). The Gaussian probability distributions for random fields Vi{x), {i = 
1,2) are given by 

P[Vi\ = jDVexp [-^lv,'ix)£ix)], (t = 0,1). (79) 

Comparing Eqs.([78|) and (^T]) , we observe the important difference between the d-wave 
superconductor and OAF. In both cases we have a continuum description in terms of 2D 
"relativistic" fermions. However, nonmagnetic impurities introduce different types of ran- 
domness to the corresponding Dirac Hamiltonians: random gauge field in the former case, 
and random mass and chemical potential in the latter case. This is related to different 
symmetries of the condensates in the two cases. 

Averaging over disorder results in the following quantum Hamiltonian of ID interacting 
fermions: 

H = J dx[r]+{-id^T3 + ^„r2)r/a + -go{: T]^r27]a -f - -gi{. vlnVa -Y]- (80) 



As before, a = 1,2, ...,r is the replica index. 

At c<j„ = 0, the case Qq = —gi would correspond to the massless chiral Gross-Neveu 
model, whose symmetry group is U{r) = U{1) x SU{r). However, since both coupling 
constants are positive by definition, this case is never realized; so even at a;„ = the model 
(pop has the 0(2r) = Z2 x S0{2r) symmetry. The important difference between these two 
groups is that the latter does not contain the continuous subgroup U{1) which survives 
in the replica limit. Therefore it is possible that at r = the model (^) has no gap less 
excitations. Thus the symmetry can be spontaneously broken, and contrary to the case of 
the d-wave superconductor, the DOS at a;„ = 0, being the corresponding order parameter, 
iVaVa), can be nonzercO. 

The breakdown of the discrete symmetry must be signalled by the development of strong- 
coupling regime in the zero-replica limit. To identify the effective couplings that are subject 
to renormalization, we rewrite the interaction terms in flSOD as follows: 



Hint = ~JJ dxi'^R,a'^R,b^aba'b''ipL,b''ipL,a' + H.C.) 

-- j dx{^p-^^^^t,b^aba'b'^L,b'tpR,a' + H.C.). (81) 

Here 

^aba'b' = ^{^aa'hw — ^ab'ha'), 

'^aba'b' = ^2Saa'Sbb' — ^l^ab'Sba' (82) 

are amplitudes of the Umklapp and momentum (chirality) conserving scattering processes, 
respectively. The bare values of the parameters in Eqs. (|8^) are 

f{0) = -{go + gi), T2{0) = -{go-gi), T,{0) = 0. (83) 

In the leading logarithmic approximation (one-loop accuracy), the renormalization group 
equations are 

§ = -[(r-i)r,-2r2]f, 
f4(f^-r?), (84) 

where / = ln(A/|e|) is a logarithmic variable. 

It is instructive to consider first a somewhat artificial case, when only backward impurity 
scattering, corresponding to a random Dirac mass m{x) in Eq.(|78|), is present, go = 0, gi ^ 0. 
The solution of Eqs. (|8^ 

1 - (r - l)gil 

indicates that a strong-coupling regime, taking place at all r > 1, is changed by a weak- 
coupling ("zero-charge") infrared behavior in the replica limit: 
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Therefore, in this case the perturbation caused by impurities is marginally irrelevant, the 
discrete 7^-symmetry remains unbroken, and the DOS will show a power-law energy depen- 
dence. 

On the other hand, when only forward scattering is considered (random chemical poten- 
tial A*(a;)), gi = 0, go ^ 0, the solution of Eqs.(84) is 

r. = o. r.(0 = -f (0 = ^^^^^ (87) 

and the situation is inverted. The weak-coupling regime, taking place at r > 1, is changed 
by a strong-coupling one at r — > 0, signalling breakdown of the discrete 7^-symmetry. The 
position of the pole of the amplitude Ti{l;r -^ 0) determines the correlation length of the 
system in the zero-energy limit, ^~^ ~ Aexp{—l/go). As a result, the DOS at c<j = is 
finite. 

It can be shown that, when both scattering processes are present, the solution of renor- 
malization group equations (^Ij) is always singular. This clarifies the special role of random 
chemical potential which drives the system towards strong coupling regime with a finite 
DOS at arbitrarily small gQ. Therefore, as long as go 7^ 0, taking into account the scattering 
between the nodes 1 and 2 would not qualitatively change this conclusion. 



VIII. CONCLUSIONS 

In this paper we have presented a detailed analysis of several two-dimensional models 
with quenched randomness. The obtained results are essentially non-perturbative and, apart 
from their practical significance, give us a new understanding of the theoretical side of the 
problem. 

We have shown that, for weak nonmagnetic impurities in a 2D d-wave superconductor, 
the whole diagrammatic series for the single-electron self energy has to be summed up. By 
explicit calculations and symmetry arguments, we have demonstrated that the DOS p(a;), 
averaged over disorder, preserves the property of being zero at zero energy. A new power-law 
behavior p(u;) ~ |u;|" was found, with critical exponent a depending on the type of included 
scattering processes. 

The feature most essential for our solution is that our models have a relativistic spectrum. 
This fact allows us to use a powerful machinery of Abelian and non-Abelian bosonization 
and the Bethe ansatz. Some of the models turned out to be exactly solvable and the main 
technical problem was to derive the replica limit of this exact solution. On the basis of 
our experience we can formulate a principle that the right procedure is to deal not with the 
energy eigenvalues, but with correlation functions. Energy eigenvalues in different symmetry 
sectors are singular at r — *^ 0, but these singularities cancel in correlation functions. 
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APPENDIX: 



Let us consider the thermodynamic Bethe ansatz equations derived for the massive sector 
of the model (|63|) in |19|. For the sake of simplicity we consider only the case N = 2. The 
free energy is described by the following system of nonlinear integral equations: 



Tp Or 

— ^ = — — / (ixcosh(7rx/2) ln{l + exp[l3er{x)]}, 

(3en = s* ln{l + exp[/3en_i(x)]}{l + exp[/3e„+i(x)]} 

—PQo COsh.{7lx/2)6n,r, 



(Al) 
(A2) 



where n 



1,2, 



and 



s* f{x) 



dy- 



f{y) 



4cosh[7r(x — y)/2\ 
In the replica limit we need to show that the following quantity has a finite value: 



lim — . 

r— ►0 T 



(A3) 



It what follows we shall consider two limiting cases: l3Qo <^ 1 and (3Qo ^1. In the first 
case, the main contribution to the free energy comes from large |x| where one can neglect 
the spectral gap and make the following approximation: 



Qo cosh(7ra;/2) ^ exp(— 7r|a;|/2 + InQo) 



(A4) 



The free energy at such /5 is a free energy of a theory with linear spectrum. It has been 
shown that this theory is the Wess-Zumino-Witten model on the SU{2) group. The free 



energy was calculated in Refs. ^ and |2J. : 



The replica limit exists: 



Fr 


TT 3r 


L 


6r + 2 " 


F 
L 


- -1 ^-' 



/5" 



(A5) 



(A6) 
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In the second case, we shall expand Fr order by order in exp(— /3Qo) and calculate the 
replica limit of this expansion. The existence of the replica limit at PQq -C 1 guarantees a 
convergence of this expansion. In the first approximation at (3Qq ^ 1 one can substitute 
€r = — cxD in Eqs. (|A2|) . The solution for n^ r has the following forma: 



l + exp[/?e„(x)] 



(n — r + 1)^ (n > r) 

(ii±i)\2 (A7) 

^\ {n<r). ^ ' 



sm - 
sin - 



Substituting e^^i into the equation for e^ we get the second iteration: 

/3e(2) = -/?Qo cosh(vrx/2) + In [2 ^^^( W^ + 2) 

sm(7r/r + 2) 

Substituting it into Eq.( |XI|) and taking the replica limit, we get the first iteration for the 
free energy: 

—— = ~P~^—r- / c?xcosh(7rx/2) exp[— (5o/9cosh(7rx/2)] 

1 r'^ I 

= --/ dpexp[-P^p' + Ql]. (A9) 

The first iteration is well defined and contains the scale Qq. It is possible to calculate the 
next term in our expansion. To do this we solve Eqs.( |A2|) linearizing them in 

dn{x) = ln{l + exp[/5e„(x)]} - ln{l + exp[/3ef ]}. 

A similar procedure was applied in Ref. |2^ for the 0(3) nonlinear sigma model. Taking the 
corresponding solutions from Appendix of Ref. ^ and substituting r = we get the third 
iteration for e,- at r — > 0: 



i/^^; 



pe'^\x) = -/3QoCosh(7rx/2) + ln(7rr) 

^JTT r-r,eM-PQocosh{ny/2)]. (AlO) 

cosh [7r(x — y)/4:\ 

From here we get the following expression for the second iteration of the free energy: 

^ = 13''^ [ dxdyexp[-Q, cosh(7ra:/2) - Q, cosh(7ry/2)] ""^|^^|"^ + y\{^ (All) 
L 8 J cosh IT [x — y)/4: 

Thus we see that the procedure does work and the r = limit exists. 
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FIGURES 

FIG. 1. Brillouin zone with the four nodes in the quasiparticle spectrum of the 2D d-wave 
superconductor labeled by 1, 1, 2, 2. Typical lines of constant quasiparticle energy are indicated. 

FIG. 2. The standard diagrammatic expansion of the self-energy T,{k,ujn) for weak static 
non-magnetic impurity scattering. The full line stands for the propagation of a quasiparticle 
with a certain momentum and each broken semi-circle corresponds to the Born scattering at the 
same impurity center. The frequency cOn is conserved and therefore omitted. 

FIG. 3. The two nodal points for the 2D orbital antiferromagnet (OAF). The Fermi surface of 
the underlying half-filled tight-binding band is the square ABCD which is also the true Brillouin 
zone for the OAF. Therefore there are only two independent nodes 1,2 in the OAF. 
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